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We consider Z2 spin liquids on the kagome lattice on the verge of a valence bond solid (VBS) 
transition, where vortex excitations carrying Z2 magnetic flux - so-called visons - condense. We show 
that these vison excitations can couple directly to the external electromagnetic field, even though 
they carry neither spin nor charge. This is possible via a magneto-elastic coupling mechanism 
recently identified by Potter et al. For the case of transitions to a 36-site unit cell VBS state the 
corresponding finite ac-conductivity has a specific power law frequency dependence, which is related 
to the crossover exponent of the quantum critical point. The visons' contribution to the optical 
conductivity at transitions to VBS states with a 12-site unit cell vanishes, however. 



I. INTRODUCTION 

Spin liquids, 2,3 albeit being Mott insulators, exhibit 
a finite ac conductivity below the Mott gap provided 
they have gapless excitations which couple to the ex- 
ternal electromagnetic field. 4-7 Recent measurements on 
the kagome material ZnCu3(OH)gCl2 (Herbertsmithite), 
which is a strong contender for exhibiting a spin liquid 
ground-state, indeed showed a specific power law fre- 
quency dependence of the ac conductivity below the Mott 
gap. 8 The nature of the ground-state in this material re- 
mains an open issue, however. While experiments are in 
favor of an almost gapless spin-liquid, 9,10 numerical re- 
sults of the antiferromagnetic Heisenberg model on the 
kagome lattice are still controversial. Projected fermion 
wavef unction studies suggest a U(l) Dirac spin liquid 
as ground-state 11 , whereas density matrix renormaliza- 
tion group (DMRG) approaches provide substantial ev- 
idence for a gapped Z 2 spin- liquid. 12-14 ; projected bo- 
son wavefunction studies also support a Z 2 spin liquid 15 . 
Large scale exact diagonalization studies show spin gaps 
and ground state energies in agreement with the DMRG 
results 16 . 

In a recent paper Potter et al. 1 identified three mech- 
anisms which give rise to a finite optical conductivity 
o{uj) of gapless spin liquids on the kagome lattice below 
the Mott gap, all leading to a characteristic a(u) ~ uj 2 
frequency dependence. They argued that a magneto- 
elastic coupling, where an applied electric field distorts 
the lattice and thereby modulates the magnetic exchange 
couplings, directly couples the external- to the emergent 
gauge field of a U(l) Dirac spin liquid and gives the 
largest contribution to the optical conductivity out of the 
three mechanisms they found. Following their approach 
we show that an external field can also couple to vison ex- 
citations of a Z 2 spin liquid by the same magneto-elastic 
coupling mechanism. Indeed, since the visons are Z 2 - 
vortices living on the dual lattice, their hopping ampli- 
tudes are modulated in accordance with the field-induced 
distortion of the direct lattice. This mechanism provides 
a direct coupling between the external electromagnetic 
field and vison excitations of a Z 2 spin liquid, which carry 



neither spin nor charge. 



In the following, motivated by the above mentioned 
DMRG studies, we focus on the situation of a Z 2 spin 
liquid 1 ' close to a valence bond solid (VBS) transition, 
where the vison gap vanishes. In a previous work 18 we de- 
rived low energy field theories for different Z 2 spin liquid 
to VBS transitions on the kagome lattice. Based on these 
results we show that the magneto-elastic coupling of vi- 
sons to the external field leads to a power law frequency 
dependence of the conductivity for transitions to a VBS 
state with a 36-site unit cell. This power law comes with 
a considerably smaller exponent than the ~ uj 2 behavior 
found in Ref. 1. By contrast, we show that the contribu- 
tion of visons to the ac-conductivity at transitions to a 
VBS state with a 12-site unit cell vanishes on symmetry 
grounds. 



Quite generally, a uniform external electric field can 
only couple to a time-reversal even spin singlet operator 
which transforms as a vector under the lattice symmetry 
group. Physically, this vector operator represents the lat- 
tice polarization P, which couples linearly to the external 
electric field E: 

5H~p.E (1.1) 

Since visons transform projectively under lattice sym- 
metry operations, we have to find an operator in terms 
of the vison fields which transforms as a vector under 
the corresponding projective symmetry group (PSG). 19 
In our previous paper 18 we showed explicitly how the 
visons transform under the PSG and we are going to 
utilize these results in order to construct such a vector 
operator on symmetry grounds. Moreover, we show that 
this operator is indeed the unique polarization operator 
which couples to the external field by the above men- 
tioned magneto-elastic mechanism. 




where the product is over an elementary plaquette of the 
dice lattice. A specific gauge choice for the Jij's which 
satisfies this frustration condition is shown in Fig. 1. 
Note that any gauge choice requires an extended unit 
cell. In Eq. (2.2) we didn't explicitly include higher order 
terms in the fields </> that are allowed by symmetry and 
describe interactions between the visons. Furthermore, 
note that the hopping amplitudes Jij are not equal to 
super-exchange amplitudes in the Heisenberg model 
above, but they are expected to be of the same order, 
since the magnetic super-exchange coupling is the only 
energy scale in the problem. 



FIG. 1: (Color online) kagome lattice (black solid lines) and 
its dual dice lattice (blue dashed lines). Shown is one ex- 
tended 12-site unit cell of the dice lattice together with a 
gauge choice of frustrated bonds (Jij — — J shown as solid 
red lines; other dice lattice bonds have Jij = J) for the ef- 
fective frustrated Ising model describing the vison excitations 
(see text). 



Important for our considerations is the fact that a dis- 
tortion of the kagome lattice inevitably leads to a dis- 
tortion of the dual dice lattice. This implies that the 
hopping amplitudes Jy of the visons on the dice lattice 
are modulated by the magneto-elastic coupling in a way 
similar to the modulation of the exchange amplitudes in 
the Heisenberg model, giving rise to a perturbation of 
the form 



II. MAGNETO-ELASTIC COUPLING 
BETWEEN VISONS AND THE EXTERNAL 
ELECTROMAGNETIC FIELD 



The basic idea behind the magneto-elastic coupling 
mechanism of Rcf. 1 is the fact that positively charged 
copper ions can be displaced within the unit cell by an ap- 
plied uniform external electric field, which in turn leads 
to a modulation of the super-exchange amplitudes be- 
tween neighboring spins. The corresponding perturba- 
tion to the Heisenberg model takes the form 



5H = 5 J^ Sj • Sj 



(2.1) 



where Si describes a spin-1/2 operator on lattice site 
i. The pattern of modulated super-exchange amplitudes 
5 Jij which couples to the external field can be inferred 
on symmetry grounds, as shown explicitly in Ref. 1. 

For our purpose, the important low energy excitations 
of a Z2 spin liquid close to a VBS transition are vortices 
carrying Z 2 magnetic flux, so-called visons. These are 
described by an effective fully frustrated transverse field 
Ising model on the dual dice lattice (see Fig. 1). Its low 
energy properties are captured by a field-theory of soft- 
spin modes 18 with a Lagrangian of the form (r denotes 
the imaginary time) 

i<j 

together with the full frustration condition 



II si S n ( J y) = - 1 

plaq. 



(2.3) 



i<3 t 



E 



(2.4) 



where P(£) is the polarization of an elementary unit cell 
of the dice lattice and the sum on £ runs over unit cells. 
We have Sij > (5ij < 0) for a squeezed (stretched) 
bond, where the hopping amplitude increases (decreases) . 
The absolute value of Sij , which parametrizes the change 
of the vison hopping amplitude with respect to the ap- 
plied electric field E can be simply estimated as 



eE 
K Cu a 



(2.5) 



with -Kcu, e, and a as effective spring constant of the cop- 
per ions, elementary charge, and lattice constant, respec- 
tively, and E = |E| is the external electric field. The pat- 
tern of modulated hopping amplitudes can be inferred on 
symmetry grounds. Following Ref. 1 we construct an irre- 
ducible representation of the modulated bonds SijJij4>i<\>j 
which transforms as a vector under the symmetry group 
of the dice lattice. 20 This operator is proportional to the 
lattice polarization operator P(£) = (P x , P y ) which takes 
the form 



Px = 



p„ 



ea 
~J 

ea 
~J 




V3 




(2.6) 
(2.7) 



where black (white) ellipses correspond to squeezed 



bonds +5Ji. 



(stretched bonds —8Jij4>i4>j) within one 



elementary unit cell of the dice lattice and we have de- 
fined 



J 



(2- 



2 





with J denoting the absolute value of the nearest neigh- 
bor vison hopping amplitude. The dimensionless param- 
eter 8 represents the magneto-elastic coupling strength 
and is typically much smaller than one. While the elas- 
tic energy KcuO 2 is on the order of leV, J takes values 
on the order of 10 meV. On the extended 12-cite unit cell 
shown in Fig. 1 the basis elements take the explicit form 



(2.9) 



= w[0 2 (R)0 5 (R)-<MR)<MR) 

+0 3 (R)0 4 (R + 2v - 2u) - fo(R)0 lo (R) 
-0 6 (R)0 9 (R) + 7 (R)0 8 (R + 2v - 2u) 
+0 9 (R)0n(R) - 0i 2 (R)02(R + 2v)] 

and similar expressions for the two other basis elements. 
Note that the sign of the gauge choice from Fig. 1 has 
been taken into account in the equation above. R = 
2m u + 2n v denotes the lattice vector of the extended 
12-site unit cell, with m,n € Z and u = (3/2, v3/2), 
v = (3/2,— v / 3/2) as the two basis vectors of the dice 
lattice. 

Also note that the specific hopping amplitudes beyond 
nearest neighbors, which had to be introduced in Ref. 18 
in order to give a dispersion to the otherwise flat vison 
bands, are not modulated here, because the distance be- 
tween these sites doesn't change. 



III. OPTICAL CONDUCTIVITY 
A. Transition to a 12-site VBS state 

We start by briefly reviewing the necessary results of 
Ref. 18. The field theory for visons at the transition from 
a Z2 spin liquid to a VBS state with a 12-site unit cell 
is an 0(4) theory with additional 4th order terms which 
break the symmetry down to GL(2,Zs). Its Lagrangian 
takes the form 



C = ((W„) 2 + (d^n) 2 +r^ + u^) 

n=l..A 

+ a ^2 ^nV'm + b [^1(^2^3 - "0204 + -03^4 



1.2 1 



+02 (01^3 



01-04 - 030'4) 
0>lV>4 + V'20'4) 
1pl(lpllh + 01^3 + ifafo)] ■ 



(3.1) 



The fields tpi are linear combinations of the soft-spin 
fields 4>i and correspond to the modes which become crit- 
ical at the VBS transition, i.e. 



,(R) 



{n) 
3 



(3.2) 



where Vj denotes the four eigenvectors corresponding to 
the highest, four-fold degenerate eigenvalue of the Fourier 
transform of the hopping matrix Jq=o a ^ momentum q = 
and j = 1, . . . , 12 is a sublattice index in the extended 
12-site unit cell (see Fig. 1 and Ref. 18 for details). Note 
that in this case 0j(R) does not depend on the lattice 
vector R of the extended 12-site unit cell. 

The four fields (?/>i, -02, 0>3, i)ja) transform under lattice 
symmetries via a four dimensional representation of the 
PSG generated by lattice translations T u , six- fold rota- 
tions Rq and reflections about the x-axis I x , which take 
the form 
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(3.3) 



We are looking for an operator which transforms as a 
vector under (3.3). Such an operator which couples to 
lattice distortions (i.e. terms in the vison Hamiltonian 
of the form 5Jij(f>i(j)j) can only be bilinear in the fields 
-0. Moreover, such an operator has to change its sign 
under rotations by 180 degrees if it is supposed to trans- 
form like a vector. Now one can readily show from (3.3) 
that (i?6) 3 = T i-e. the fields transform back to them- 
selves after 180 degree rotations, thus it is not possible 
to construct a vector operator out of the fields ipi with- 
out invoking spatial gradients. Since the external electric 
field is homogeneous and only couples to the zero mo- 
mentum component of such a vector operator, any PSG 
invariant bilinear term in the 0j's involving gradients is 
ruled out. In principle it would be possible to construct 
a fourth order term involving gradients which transforms 
as a vector and gives a non-zero contribution to the con- 
ductivity when coupled to a homogeneous external field, 
but such an operator cannot correspond to the lattice 
polarization operator as it is not a bilinear. Since we 
are not aware of a mechanism which couples the electric 
field to a quadrilinear operator in the fields ip, we do not 
pursue this route further at the moment. 

We can also directly show that the lattice polarization 
operator vanishes for the case of a transition to a 12- 
site VBS state by using Eqs. (2.6), (2.7), (2.9) and (3.2). 
Since 0j(R) doesn't depend on the lattice vector R we 
get 
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where the last equation follows after using the explicit 
form of the eigenvectors Vj from Ref. 18. Similarly also 
the other two basis elements that are used to construct 
the polarization operator shown in Eqs. (2.6) and (2.7) 
vanish. Consequently it is not possible to construct a lat- 
tice polarization operator which couples to the external 
electric field via the magneto-elastic mechanism. In fact 
we see that the bonds related by a 180 degree rotation 
cancel each other, confirming our earlier argument. 



where q = ±Qi = (0, are the momenta where 

the vison dispersion has its minima. When inserting Eq. 
(3.14) into Eqs. (2.6) and (2.7) it is important to keep 
only the zero-momentum components (i.e. the terms in- 
dependent of the lattice vector R) since the external elec- 
tric field is homogeneous and will only couple to such 
terms. Using the shorthand notation 



E 

n=1...4 



(3.15) 



we obtain 



B. Transition to a 36-site VBS state 

For a transition to a VBS state with a 36-site unit 
cell the visons are described by an 0(8) theory with ad- 
ditional fourth and sixth order terms which break the 
symmetry down to GL(2,7L-$) x D3. The explicit form of 
the Lagrangian and the corresponding eight-dimensional 
matrix representation of the PSG can be found in Ref. 18. 
In contrast to the previous case, we can straightforwardly 
construct a homogeneous bilinear operator which trans- 
forms like a vector under the PSG by making the ansatz 

(r = (v x ,p y )) 



P, = > Oiji/jit/jj, V y = hjipiipj (3.5) 

f,3=1...8 



E 

i,j=1...8 



and determining the coefficients of the matrices (a) 



ij — b{j via 
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(3.11) 



This vector operator is uniquely determined up to a con- 
stant prefactor and takes the form 



(3.12) 



(3.13) 



As we will show now, this operator corresponds to 
the unique lattice polarization operator P = ea5 T in 
Eqs. (2.6) and (2.7). For transitions to a VBS state with 
a 36-site unit cell, the soft-spin modes <fi are related to 
the critical modes ip via 18 
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(3.14) 
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The exponential factors arise from bonds connecting sites 
in adjacent unit cells. Using this and similar expres- 
sions for the two other basis elements, we reproduce Eqs. 
(3.12) and (3.13) exactly up to a multiplicative constant. 

The frequency dependence of the conductivity can now 
be obtained using Kubo's formula 



(P.P 



(3.17) 



where the frequency dependence of the polarization cor- 
relation function is determined by the scaling dimension 
of the polarization operator P^. In Ref. 18 we couldn't 
find a stable fixed point at one-loop order for our 0(8) 
theory including the additional fourth order terms. Mak- 
ing the simplifying assumption that the critical point be- 
tween the Z2 spin-liquid and the VBS phase is described 
by the standard 0(8) Wilson-Fisher fixed point, the po- 
larization operator in Eqs. (3.12) and (3.13) can be ex- 
panded in terms of components of the traceless symmet- 
ric bilinear tensor operator of the 0(8) model 



T mn (x) = ip n (x)ip m (x) 



(3.18) 



where ip n (n £ {1, . . . , 8}) now represents the real- and 
imaginary parts of the four complex fields in Eqs. (3.12) 
and (3.13). We denote the scaling dimensions of this 
tensor operator, and thus the polarization operator, by 
At- Then the frequency scaling of the ac conductivity 
in d = 2 dimensions is given by 



,2A T -2 



(3.19) 



The scaling dimension can be expressed in terms of the 
crossover exponent 4>t by Ay = 3 — 4>tI u i and has been 
evaluated to six loops 21 as 4>t = 1.55. Using v = 0.830 22 



4 



we obtain 

2A T - 2 = 0.27 (3.20) 

This exponent is considerably smaller than two, giving 
rise to a large optical response compared to the mecha- 
nisms discussed in Ref. 1 at frequencies below the mag- 
netic super-exchange coupling. The resulting expression 
for the frequency dependent conductivity finally takes the 
form 




(3.21) 



where we have reintroduced Planck's constant h. Note 
again that J denotes the nearest neighbor vison hopping 
amplitude here, which is expected to be on the order of 
the magnetic super-exchange coupling of the underlying 
Heisenberg model, since this is the only energy scale in 
the problem. 

IV. DISCUSSION 

We calculated the optical conductivity of Z 2 spin liq- 
uids at a quantum critical point to a VBS state with a 
36-site unit cell, and found a power-law frequency depen- 
dence of the conductivity ct(cj) ~ w 27 , a particularly 
small exponent. This indicates a large optical conduc- 



tivity at low frequencies at the critical point, and the 
enhancement should also persist away from the critical 
point. 

This calculation is based on the simplifying assump- 
tion that the transition is described by the Wilson-Fisher 
fixed point of the 0(8) model. It is still an open question, 
if the full theory including the 0(8) symmetry break- 
ing terms presented in Ref. 18 exhibits a stable fixed 
point, and it would be interesting to study its critical 
behavior in the future. Going beyond a one-loop calcu- 
lation is a daunting task, however, given the complicated 
form of the Lagrangian (see Eq. (3.24) in Ref. 18). A 
viable approach would be to study the frustrated Ising 
model directly using Monte-Carlo simulations. In any 
case, the smallness of the optical conductivity exponent 
relies mainly on the fact that polarization operator is a 
bilinear in the field of the 0(8) model, with no addi- 
tional spatial gradients, and so we can reasonably expect 
a small exponent at any possible fixed point which breaks 
the 0(8) symmetry down to GL(2,Z 3 ) x D 3 . 
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